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Abstract. We review here some universal aspects of the physics of two-electron molecular transistors in the absence of strong spin-orbit 
effects. Several recent quantum dots experiments have shown that an electrostatic backgate could be used to control the energy dispersion 
of magnetic levels. We discuss how the generically asymmetric coupling of the metallic contacts to two different molecular orbitals can 
indeed lead to a gate-tunable Hund's rule in the presence of singlet and triplet states in the quantum dot. For gate voltages such that 
the singlet constitutes the (non-magnetic) ground state, one generally observes a suppression of low voltage transport, which can yet be 
restored in the form of enhanced cotunneling features at finite bias. More interestingly, when the gate voltage is controlled to obtain the 
triplet configuration, spin S = 1 Kondo anomalies appear at zero-bias, with non-Fermi liquid features related to the underscreening of a spin 
larger than 1/2. Finally, the small bare singlet-triplet splitting in our device allows to fine-tune with the gate between these two magnetic 
configurations, leading to an unscreening quantum phase transition. This transition occurs between the non-magnetic singlet phase, where 
a two-stage Kondo effect occurs, and the triplet phase, where the partially compensated (underscreened) moment is akin to a magnetically 
"ordered" state. These observations are put theoretically into a consistent global picture by using new Numerical Renormalization Group 
simulations, taylored to capture sharp finie-voltage cotunneling features within the Coulomb diamonds, together with complementary out- 
of-equilibrium diagrammatic calculations on the two-orbital Anderson model. This work should shed further light on the complicated puzzle 
still raised by multi-orbital extensions of the classic Kondo problem. 
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1. Introduction 

Quantum dots, artificial nanostructures with quantized 
electronic charge that can be probed by electrical 
transport, offer promising ways of manipulating the spin 
of electrons in atomic size devices [H [2]. By using 
a combination of lithographic methods, electrical gates 
and applied magnetic field, a high degree of control of 
single-electron magnetism can be thus achieved, allowing 
e.g. the realization of spin-qubits [3 . By opening the 
tunneling barrier between quantum dots and the metallic 
electrodes used as contacts, quite fascinating physics 
emerges at strong tunnel coupling. In the case of an 
odd number of trapped electrons in a dot with a single 
orbital well separated in energy from other excitations (due 
to confinement and Coulomb blockade), an artificial and 
tunable version of the spin 5* = 1/2 Kondo effect can be 
realized [U |5j |6]. In that situation, magnetic screening 
of the single spin by the Fermi sea occurs below the so- 
called Kondo temperature T^, giving rise to a conductance 



increase up to the maximum unitary value 2e^//i, with e 
the electron charge and h Planck's constant [7l|8]. The first 
observation in semiconducting quantum dots of the Kondo 
effect [9l[T0] has triggered intense experimental research to 
make similar observations in other physical systems such 
as carbon nanotubes [11 and molecular devices dU [13], 
see Ref. [14 for a review. In parallel, many theoretical 
works have followed, and at present reliable methods, 
such as the Numerical Renormalization Group (NRG) 
[151 US? have been developped. These calculations can 
allow quantitative comparison to experiments, especially 
because the spin S = 1/2 Kondo effect shows simple 
universal scaling laws, that can be checked in quantum 
transport measurements. At present, one important goal 
is to understand in similar detail the electronic transport 
through multi-orbital dots, a situation clearly raised both 
by semiconducting and molecular devices. The next step 
towards complexity, two-electron quantum dots, seems 
deceptively small, yet brings a richness that is still not 
fully clarified. Theoretically, a multitude of new effects 
have been predicted over the years in the realm of Kondo 
physics, and many have been observed thanks to the 
great tunability of quantum dot systems. We briefly 
give an overall review of this physics in the rest of this 
introduction, with additional details provided in Section |4] 

Let us first focus on the case where the two orbitals 
experience strong ferromagnetic Hund's rule, so that the 
triplet configuration constitutes the main low lying states. 
This can result in the spin S = 1 underscreened Kondo 
effect if a single screening channel is active, as initially 
proposed by Nozieres and Blandin [17], and only recently 
experimentally observed ^8] [19] . This rather exotic Kondo 
effect shows a singular (logarithmic) approach to strong 
coupling, in contrast to the regular Fermi liquid behavior 
of fully screened Kondo impurities observed in odd-charge 
spin S = 1/2 quantum dots Because Hund's rule 

competes with the orbital level spacing, an intra-orbital 
singlet state is usually close (and often lower) in energy to 
the triplet states, and the competition between singlet and 
triplet can give rise to very rich physics. 

One interesting example is the so-called singlet- 
triplet "transition" , which realizes a fully screened Kondo 
effect by bringing into degeneracy singlet and triplet 
magnetic levels (it is thus rather a crossover than a real 
phase transition), for which two different scenario can 
be obtained. The simplest situation occurs in carbon 
nanotubes (or molecular quantum dots), where a strong 
Zeeman effect is used to cross the singlet state with the 
lowest triplet state in the presence of a single screening 
channel |20l [21J , showing the emergence of a strong Kondo 
enhanced conductance at the crossing point. Alternatively, 
similar physics can be realized in semiconducting vertical 
quantum dots, where the magnetic field can be used to tune 
by orbital effects the splitting between the singlet and the 
three-fold degenerate triplet [22] [23l [24l [25] (the Zeeman 
effect can be neglected due to the small g-factor). Here the 
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vertical structure preserves the orbital quantum numbers 
during the tunneling process, so that two screening channel 
are active, and the high spin states are again fully Kondo 
screened at the crossing point. 

One can then naturally wonder about the fate of 
a singlet to degenerate triplet crossing in the presence 
of a single screening channel [26} l27t 128} [29l [3Ql l3T| 
[32| [33| [34] , which can in principle be expected for 
lateral quantum dots [35l |36^ and molecular devices 
[34 . Interestingly, this situation gives rise to both 
underscreening physics [171 HE] and an unscreening singlet- 
triplet quantum phase transition (QPT). In that case, a 
zero-temperature phase transition can take place between 
a non-magnetic singlet state and a partially compensated 
spin S = 1 (giving an remanent spin 5 = 1/2 in the ground 
state). The transition is really sharp at zero temperature 
as the entropy changes from zero (in the singlet phase) 
to log(2) in the underscreened Kondo phase. This is 
in contrast to the so-called singlet-triplet "transition" 
induced by a magnetic field, where the entropy vanishes 
in both "phases" (as discussed in the previous paragraph). 
Several different theoretical proposals have also been made 
for the unscreening QPT ^3 EH E^, which was 
shown to be even robust to valence fiuctuations. Our 
detailed experimental observation [34] of this singlet-triplet 
quantum phase transition in a two-electron molecular 
quantum dot will be one of the main topics of this review, 
and the various scenarios presented above will be reviewed 
with more details in the concluding Section [4] 

Out-of-equilibrium Kondo phenomena have also been 
investigated in two-electron quantum dots, such as Kondo- 
enhanced cotunneling lines at finite bias 03 [Ml [40] , 
and these will be mentioned later on. We stress also 
that spin-orbit interaction has been given some recent 
attention both experimentally 2Tl|42][l2 and theoretically 
[43l|44l|45l[l9l|46l[47] for even-charge quantum dots, as weh 
as for odd-electron quantum dots [48 , and these effects 
could possibly play an important role for the interpretation 
of some experiments made in semi-conducting quantum 
dots [35l [36]. We also remark that carbon nanotube 
quantum dots can show a wealth of interesting Kondo 
states due to the enhanced SU(4) symmetry, both in 
single and doubly occupied situations, but these system 
go beyond the present discussion (see e.g. Ref. [49l [50] 
and references therein). 

A different physical situation occurs when the two 
orbit als interact with antiferromagnetic exchange. In the 
case of a single screening channel, a two-stage Kondo effect 
occurs for large Kondo coupling compared to the interspin 
exchange. Indeed, after complete Kondo screening of 
the first spin, the resulting Fermi liquid state is able to 
absorb the second spin in a second stage of screening 
[51]. The situation of two screening channels leads 
potentially to even more exotic physics. This provides a 
realization of the so-called two-impurity Kondo problem 
[52] , which realizes a competition between Kondo screening 



and (RKKY) exchange and presents a non-trivial quantum 
phase transition. The experimental observation of this 
transition turns out to be difficult f53", "ST despite previous 
attempts [55], due to relevant physical processes at the 
non Fermi liquid fixed point (essentially the interdot 
tunneling). 

Finally, it is also worth noting the complicated role 
of charge when capacitive coupling between two dots (or 
two orbitals) is relevant [56^, which can give rise to complex 
charge skipping in the filling scheme of the dots, a situation 
not fully elucidated in the presence of spinful electrons. 

On a mathematical level, all the important complexity 
of these various phenomena can be expressed in complete 
generality by a two-orbital Anderson impurity model, 
which contains a great deal of physical parameters 
compared to the single orbital Anderson or Kondo 
Hamiltonians. One crucial question is thus to understand 
what classes of universal behaviors can be expected in such 
a complicated situation. Going one step further, we also 
mention that three-orbital models show even more complex 
and interesting physics, such as non Fermi liquid fixed 
points, a topic that goes far beyond the present review 
(see Ref. [53 |58l [59l EQI [SU for a few recent theoretical 
studies and Ref. [62] for a review). 

On an experimental point of view, multi-orbital 
effects are always present, but their relevance will depend 
on many factors, most primarily the energy spacing 
between the two relevant orbitals, and the structure of 
the hybridization matrix. It is clear that disentangling 
the mechanisms at play from transport measurements 
only is not an easy task, yet many useful information 
can be gathered from finite bias cotunneling spectroscopy, 
possibly in the presence of magnetic fields. The present 
study will illustrate this general strategy, and show how 
the magnetic states involved can be identified through an 
analysis of Zeeman split cotunneling lines. 

The paper is organized as follows. Section [2] will 
present recent observations made with fullerene quantum 
dots obtained by electromigration. These will concern 
finite bias cotunneling spectroscopy of magnetic states, 
gate-dependent tuning of singlet and triplet states, the 
spin S = 1 underscreened Kondo effect, and the 
unscreening quantum phase transition. Connection to 
related results in both molecular and semiconducting 
nanoelectronics will also be made, showing the wide 
range of application of those ideas. Section [3] will 
discuss modelization of the problem, and the combined 
use of near-equilibrium Numerical Renormalization Group 
calculation and diagrammatic methods in the non- 
equilibrium situation to reproduce semi-quantitatively 
most of the features seen in the experiment. Finally 
section [4] will close the paper by presenting a more detailed 
review of the various interesting Kondo phenomena 
involving singlet and triplet excitations in two-electron 
quantum dots, that will set our observations in a broader 
perspective. 
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2. Experiments in even charge quantum dots: 
gate control of spin states and underscreened 
Kondo effect 

Kondo physics has been widely reported in a variety 
of quantum dot systems, based on semiconducting 
heterostructures [9l [10], carbon nanotubes [11] and 
molecular devices [HI [131 (311 [HI [H] , showing the great 
universality of this phenomenon. Usually the Kondo effect 
is best observed for odd charge quantum dots, because 
the formation of a net spin S = 1/2 is always guaranteed 
in that case (for even charge dots, the ground state 
turns out often to be a non-magnetic singlet). While 
Kondo signatures associated to a triplet configuration in 
two-electron quantum dots was previously observed [36] . 
little attention was devoted to the detailed study of this 
phenomena. One reason is that the Kondo temperature for 
spin = 1 is generically quite low, making semiconducting 
systems not suitable to explore this physics. In contrast, 
molecular devices display much bigger energy scales, 
allowing to investigate in great detail the richness of the 
Kondo effect in two-electron quantum dots. The aim of 
this paper is to expose the generic and if possible universal 
features in such two-electron systems, both at the light 
of experiments on fullerene transistors and also thanks 
to state-of-the-art many-body simulations. This section 
will be devoted to the presentation and discussion of the 
experiments, while the theory will be presented in Sec. [3] 
In our previous experimental work [34 , a gate-tuned 
transition from a low zero-bias conductance state to a high 
zero-bias conductance state was reported in a molecular 
quantum dot containing two electrons, see Fig.[l] The data 
also show a striking collapse of the magnetic excitations 
on both sides of the transition, and can thus naturally 
be explained by the so-called singlet-triplet unscreening 
quantum phase transition [26^ ^27^ 28^ 22l [3Ql |32] , where the 
low conductance regime (left part of lower panel in Fig. [T]) 
corresponds to a singlet state binding of the molecular 
orbitals, while the high conductance regime (right part of 
lower panel in Fig. [T]) is associated to the underscreened 
Kondo effect of a spin 5* = 1 triplet. We note that similar 
features of Kondo anomaly collapse have been reported in 
other two-electron quantum dot systems |36l [63] [M] [19] , 
but did not show a concomitent crossing of the magnetic 
excitations. Recent theoretical work [19, l45l proposed 
spin-orbit effects (or more generally magnetic anisotropics) 
as a possible alternative interpretation for some of these 
other experiments. We will show however that spin-orbit 
coupling can be neglected in our particular experiment, 
which is a crucial requirement [T9^, '^S for the robustness 
of our observation of Kondo underscreening [18^ . Finally, 
recent experimental and theoretical studies (65] [66] [671 [181 
|68[ l69[ [70] have confirmed several microscopic mechanisms 
at play in the singlet-triplet quantum phase transition, 
such as the gate tunable Hund's rule, giving extra strength 
to our initial interpretation. We will refer the reader to our 
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Figure 1. Observation of a Kondo screening-unscreening 
transition in an even charge Coulomb diamond of the 
molecular transistor [34]. The lower panel is a close up inside 
the charge N = 2 Coulomb diamond seen on the upper panel. 
Singlet /triplet states of the molecule are determined from the 
absence/presence of a zero bias Kondo anomaly, as well as by 
the Zeeman spectroscopy analysis shown in Fig. [2] 

previous work [34l [7T| [T8] for the experimental details, 
and will present in what follows a physical discussion of 
our observations. 

2.1. Cotunneling spectroscopy: nature of the magnetic 
ground state and excitations 

Identifying the molecular magnetic states involved on 
the quantum dot is certainly a key step in the further 
interpretation of our transport experiments. Because 
Kondo physics results at relatively large tunnel coupling to 
the electrodes, sequential tunneling lines, occuring along 
the edges of the Coulomb diamonds, are typically very 
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broad and masking the useful information on the ahowed 
transitions from a charge state to the next. For this 
reason, cotunnehng spectroscopy done within the Coulomb 
diamonds (hence at fixed charge on the quantum dot) 
turns out to be a very useful tool. Because cotunnehng 
spectroscopy reveals transition lines between states with 
the same charge, their intensity is much weaker than 
transition with a valence change occuring at the edges of 
the Coulomb diamonds. This strategy has clear advantages 
for improving the spectral resolution, but requires at the 
same time careful experimental measurements, performed 
at subKelvin temperatures and using efficient filtering of 
the electric noise. 

Let us thus focus on the cotunnehng lines seen in the 
even charge Coulomb diamond of our experiment (lower 
panel of Fig. [T]), and study finite bias spectroscopy in 
a magnetic field. If a spin triplet is really the ground 
state for large positive gate voltage, as is hinted by the 
presence of a zero-bias anomaly, it should be split in 
an obvious manner by the Zeeman effect, and magnetic 
transitions should appear accordingly. The transport data 
shown on the right panel of Fig. [2] clearly confirms this 
interpretation, as the selection rules from the ground 
triplet states to the excited singlet state are precisely 
the ones expected. Similarly, magnetic-field induced 
transition from the singlet to the lowest triplet is clearly 
evidenced on the Zeeman lines obtained as a function 
of the magnetic field at fixed gate voltage (left panel of 
Fig.[2|. This understanding of which states are involved in 
the experiment provides already useful information for the 
interpretation and modelization presented in the following. 
Also, the absence of anticrossings between the magnetic 
excitations shows that singlet and triplet remain pure spin 
states. This points to the absence of spin-orbit effects 
(see Ref. ^TJ 02 where such effects are evidenced in a 
self- assembled semiconducting quantum dot and carbon 
nanotubes respectively), which is for instance a crucial 
point for the robustness [19 , 45 of an underscreened Kondo 



effect (discussed in section 2.3). 



2.2. Driving mechanism for the gate- controlled 
singlet-triplet transition 

One striking question should arise at the light of the above 
analysis: why a magnetic transition can be driven by a 
change of the electrostatic potential? Such phenomena, 
not fully elucidated in our original work has received 
a lot of experimental and theoretical attention recently 
[65l [Ml EH EHl ES] and has a simple explanation in 
terms of molecular levels renormalized by virtual charge 
fluctuations. 

A crucial piece of this puzzle was infact already given 
in an observation made in our previous work [34], namely 
that under a magnetic field the crossing of the singlet 
|0) and the lowest l) triplet does not lead to an 
enhanced conductance (as seen in the left panel of Fig. 
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Figure 2. Identification of the singlet and triplet excitations 
by the Zeeman effect in our device [34]. The left panel shows 
the magnetic field dependence of the cotunnehng lines at a 
fixed gate voltage taken in the singlet region. The right panel 
presents a gate voltage scan for a given value of the magnetic 
field (3 Tesla). In both cases, spin transitions follow the 
selection rules associated to a singlet-triplet scenario. 



|2|, in stark constrast to the (single-channel) singlet-triplet 
"transition" scenario [20*, "27. This implies that, at the 
level crossing, the degenerate "doublet" {|0), l)} 
induced by the magnetic field possesses a vanishingly low 
Kondo temperature. This necessarily means that the 
two molecular orbitals (1,2) must have an asymmetric 
tunneling amplitude to the leads, namely ti ^ ^2, since 
the effective Kondo exchange interaction associated to this 
pair of states is proportional to tit2/U (with U the dot 
charging energy). Indeed, a "spin-flip" between the singlet 
and lowest triplet configurations involve tunneling out of 
orbital 1 and tunneling into orbital 2 (and vice- versa). 

Now a qualitative argument, analogous to the one 
proposed in Refs. [65^, ^6 for odd charge quantum dots, 
and based on detailed perturbative calculations in the 
tunneling amplitude [28l [31], shows that the existence 
of this orbital tunneling asymmetry is the key to the 
gate-controlled singlet-triplet splitting. The basic idea 
relies on a lowering of the singlet (resp. triplet) energy 
on the left (resp. right) side of the N = 2 Coulomb 
diamond by virtual charge fluctuations to the = 1 
(resp. A" = 3) neighboring Coulomb diamond. This is 
due to the fact that the flrst orbital has indeed a higher 
probability of tunneling than the second one (see Fig. 
|3|, so that an energy gain of order (t^ — tD/U can be 
achieved by either the singlet or the triplet conflguration, 
depending on the proximity to the left or right Coulomb 
diamond respectively. This hypothesis is in agreement 
with our previous observation [34] of asymmetric tunneling 
amplitudes, as discussed above, and this simple mechanism 
explains how a gate-tunable Hund's rule can be achieved 
in our device. 
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Figure 3. Panel a: stabilization of the singlet state on the left side of the N = 2 diamond via virtual fluctuations to the = 1 
diamond. Panel b: stabilization of the triplet state on the right side of the N = 2 diamond via virtual fluctuations to the N = 3 
diamond. The mechanism is effective when level 1 is more strongly hybridized than level 2 to the leads. See Refs. [28l for a 
detailed theoretical treatment by perturbation theory in the tunneling. 
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Figure 4. Illustration of the gate-tuning of magnetic states and the bending of the cotunneling lines at the diamond edges. Panel 
a: our single molecule transistor data for an even charge N = 2 quantum dot [34| where a gate-tunable Hund's rule is achieved. 
Panel b: carbon nanotube device by Hauptmann et al. [66| for an odd charge N = 1 quantum dot coupled to ferromagnetic leads, 
where a gate-tunable Zeeman splitting is realized. 



On a more quantitative level, one can also understand 
why the dispersion of the magnetic states becomes non- 
linear near the edges of the Coulomb diamond, as seen 
in Fig. [4j Lowest order perturbation theory in the 
tunneling events associated to the first orbital provides a 
correction to the bare singlet-triplet splitting, that diverges 
logarithmically near the diamond edge [66]. This sharp 
enhancement is clearly responsible for the bending of the 
magnetic excitation lines near the edge of the Coulomb 
diamonds, and should be expected to be a generic effect. 
Indeed, such features have now been widely reported, 
and for illustration purposes we present data on carbone 
nanotube quantum dots [66] , taken in the odd charge sector 
for a device with ferromagnetic electrodes, see Fig. ^p. 
There, the tunneling asymmetry concerns rather the up 
and down states of the impurity spin S = 1/2 (so that 
e.g. and in a similar fashion to our findings, 

leads here to a gate-tuned Zeeman effect for these magnetic 
excitations (the non-linearity near the diamond edges is 
also quite apparent in this data set). 

One can thus infer that a gate-induced singlet- 
triplet transition can be considered a general effect that 
will however be observed according to some quantitative 
criterium: the bare singlet-triplet splitting of the isolated 
molecule (that results from the competition of level 
spacing. Coulomb repulsion and Hund's rule) must be 



smaller than the singlet to triplet energy dispersion that is 
introduced by the asymmetric coupling to the electrodes. 
The singlet-triplet transition can thus only occur at 
relatively strong hybridization to the leads. In alternative 
cases where the singlet-triplet splitting turns out to be 
larger than the hybridization shifts, the quantum dot 
keeps intact its magnetic ground state (singlet or triplet) 
throughout the whole Coulomb diamond. This is for 
instance routinely seen in less strongly coupled quantum 
dot systems, such as carbon nanotubes [39], or our own 
measurements in other Ceo devices [18]. 

2.3. Triplet side: spin S = 1 underscreened Kondo effect 

If one focuses on the side of the transition where a Kondo 
ridge is visible (large positive gate voltage), a clear spin 
S = 1 Kondo anomaly is observed. One interesting 
question is the nature of this Kondo state, which will 
crucially depend on the number of active screening 
channels (at the temperature where the experiment is 
performed). In the case of two screening channels, 
the spin 5 = 1 is fully quenched [29 , resulting in 
a conventional Fermi liquid ground state, while for a 
single screening channel, the spin 5* = 1 is partially 
compensated. This underscreening process is associated 
to logarithmic deviations to Fermi liquid theory [17]. We 
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Figure 5. Comparison of the (dimensionless) conductance 
measured in odd and even Coulomb diamonds (squares) 
to NRG calculations (lines) [18| , in units of the Kondo 
temperature Tx- 



note that Kondo anomalies with even charge states where 
previously observed, although quite rarely [45^, in some 
other semiconducting quantum dot experiments [72j [73] , 
and were not studied in great detail. 

Here, we argue that molecular quantum dots will 
generically end up in a single screening channel situation, 
in the sense that the Kondo coupling to a second screening 
bath, although always present, will in general be relatively 
small (compared to the most strongly coupled channel) , so 
that the associated Kondo temperature for a second stage 
of full screening of the impurity spin will be exponentially 
small. Several arguments are in favor of this situation: 
i) tunneling from the electrodes is typically monomode 
because the electromigration process is stopped at the 
breaking up of the metallic contacts; ii) coupling of the 
molecule to the source and drain is also typically quite 
asymmetric (conductance maxima of Kondo anomalies in 
molecules range from a few per-thousands to half the 
conductance quantum [13 ); iii) the two molecular orbitals 
involved can also be asymmetrically coupled to the leads. 



as is seen at least in our device, see discussion in section 2.2 



To explain this more presicely, let us introduce the 
matrix of tunneling amplitudes ^29^ 



tL2 tR2 



(1) 



where we describe our molecular transistor by two orbital 
levels (1,2) coupled to two metallic leads (L, R) with single 
electronic mode. A screening channel (A =even,odd) is 
associated to each eigenvalue of this matrix, from which 
antiferromagnetic Kondo couplings Jx = 4|tAp/^add 
result between the localized orbitals and the conduction 
electrons (£^add is a measure of the addition energy on 
the molecule). Typically the even combination of the 
wavefunctions results in a larger tunneling term than 
the odd one, namely teven > ^odd, because of generic 
source/drain (L/R) asymmetries, or also orbital level (1/2) 



asymmetries. In turn, the two exchange interactions will 
show relatively different magnitudes, Jeven > ^odd^ so that 
a first stage of Kondo underscreening of the 5* = 1 triplet 
will occur at temperatures T^®" cx exp(— 1/po^even), with 
po the leads density of states. A second stage of screening 
of the remanent doublet state will develop only at the 
much lower temperature T^^^ oc exp(— 1/po-^odd) < Tf^^^, 
owing to the exponential dependence of the Kondo scale 
in the exchange interaction. This implies that spin S = 1 
molecular quantum dot experiments covering a couple of 
decades in temperature will only show underscreening, 
unless fine-tuning of the tunneling matrix is realized. 

This discussion can be put onto firm ground 
by making quantitative comparisons of the measured 
temperature-dependent conductance to Numerical Renor- 
malization Group calculations, see Fig. [5] and Ref. for 
further details. The slow saturation of the low temperature 
conductance for the even charge spin = 1 regime, asso- 
ciated to the logarithmic approach to the underscreening 
fixed point, is readily contrasted to the rapid Fermi liquid- 
like behavior at low temperature for the odd charge spin 
6* = 1/2 Kondo effect. This quantitative comparison gives 
also more convincing evidence for our observation [34 of a 
quantum phase transition when singlet and triplet states 
meet at the critical gate voltage. 

2.4' Singlet-triplet unscreening quantum phase transition 

The sudden disappearance of the underscreened Kondo 
resonance in the middle of the even charge Coulomb 
diamond (see lower panel in Fig. [T]) is the sign of a 
possible quantum phase transition between two different 
sorts of ground states. From the identification of the 
singlet and triplet magnetic states performed in section[2]T] 
and the confirmation that a single screening channel 
is active (see discussion in sections 2.2 and 2.3), the 



obvious scenario to follow is the so-called singlet-triplet 
unscreening transition [26l |27l |28l El |32] . In this picture, 
the transition is driven by the merging of singlet and 
(underscreened) triplet magnetic excitations together at a 
quantum critical point. The transition is continuous, but 
shows an entropy change at the zero-temperature critical 
point: the non-degenerate singlet state cannot be smoothly 
related to a partially screened spin with log (2) remanent 
entropy. Because the zero-temperature limit where the 
transition takes really place is unaccessible to experiments, 
we discuss the associated signatures of this transition at 
finite temperature. We already note that the physics 
here is more complex than a simple singlet /triplet (or 
singlet /doublet) level crossing, because of the electronic 
correlations arising from the electrodes, and cannot be 
described by a first order transition. 

The upper panel of figure [6] shows previously unpub- 
lished results for the non-linear conductance measured for 
several gate voltages across the transition. Far on the sin- 
glet side (for gate voltages Vg <1.87V), a well formed gap 
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Figure 6. Differential conductance as a function of bias 
voltage for several gate voltages across the transition. In both 
plots, arrows denote the dispersive singlet/triplet excitations. 
The upper panel shows the experimental measurement [34| 
and the lower panel the theory (68] [69] performed in 
section |3.2| The theoretical calculations are performed for 
several values of across the QPT, with parameters D = 20, 
Et =0, Es = -0.04, Ti = 1.106, r2 = 0.38, and T = 10"^. 
In the experimental plot (upper panel), the four lower curves 
(associated to low conductance at small bias) are taken on the 
singlet side, while the upper curves correspond to the triplet 
side. 



opens up in transport, showing the strong binding of the 
molecular levels into a tight singlet unit. Yet an out-of- 
equilibrium enhancement of the cotunneling lines (to the 
triplet excitations) is witnessed [71 , similar to previous ob- 
servations made in carbon nanotubes [39] . Upon approach- 
ing the transition (for gate voltages 1.88V< Vg <1.9V), the 
finite bias shoulder progressively disperses towards zero 
bias, while a slower logarithmic-like decrease of the con- 
ductance is now observed instead of a hard singlet gap. 
These features can be explained within the singlet-triplet 
scenario [26l |27l |28l 132] by a two-stage Kondo process for 
the formation of the singlet ground state, as demonstrated 
by the complex temperature dependence of the non-linear 
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Figure 7. Differential conductance on the singlet side of 
the transition for several temperatures from 40mK to 15K, 
showing the two-stage kondo process. 



conductance, see Fig. [7| The physical interpretation of 
these data is the following ^|2i|29l[30l[3ll|3a|33[34]: 

at high-temperatures, the two-orbit als act basically as two 
decoupled spin S = 1/2, and the most strongly coupled 
spin undergoes a first screening process associated to a 
Kondo scale Tx^s=i/2i while the more weakly coupled spin 
remains essentially free. This is witnessed in Fig. [T] by the 
formation of a standard Kondo anomaly (with maximum 
conductance at zero bias) for temperatures T > 600mK. 
By further cooling, interaction between the remaining un- 
screened spin and the Fermi sea builds up via the already 
screened S = 1/2 moment, and leads to a second stage of 
screening that fully quenches the total entropy. By adia- 
batic continuity to the singlet state realized far from the 
transition, the conductance can be shown to drop to small 
values, so that a Kondo-like dip is seen in transport asso- 
ciated to a second Kondo scale T*. The presence of two 
different energy scales is clearly shown in Fig. [7| while the 
crossover from strong singlet binding to the Kondo-like sin- 
glet binding appears by the change from a hard gap to a 
softer pseudogap in Fig. [6| 

On the triplet side of the transition, also dispersive but 
less pronounced finite-bias excitations to the singlet state 
are clearly seen in Fig. [6] These features are also observed, 
albeit more mildly, in the temperature-dependent zero-bias 
conductance by a distinctive shoulder for temperatures of 
the order of the singlet excitation threshold, see upper 
panel in Fig. |8] Also quite striking in the experimental 
data on the triplet side is the ongoing increase with 
temperature of the zero bias conductance, associated to 
the underscreened Kondo effect discussed previously. We 
now turn to the modelization and the microscopic analysis 
of this physics, giving strength to the interpretation of 
the experiment in terms of the unscreening singlet-triplet 
quantum phase transition. 
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3. Theoretical analysis of the singlet-triplet 
unscreening transition 

3.1. Two-orbital Anderson model 

A quite generic model for two-electron quantum dots 
includes both relevant levels and their hybridization with 
the leads (see e.g. Ref. [38 ): 

H = i^dot + Pleads + ^mix- (2) 

The term iYdot above describes two energy levels in the 
quantum dot and their mutual interactions: 

^dot = ^ ^jriia^^ Uini^nii^Ui2nin2-JHSiS2, (3) 



where n^^ = dl^di^, rii = ^^rii^, Si= X^^/? 4a^a/3<^i/3/2, 
and dl^ creates an electron with spin a on the dot 
level i =1, 2. The first term in Eq. Q describes 
the single-particle energy of both levels, the second and 
third terms are the intralevel and interlevel Coulomb 
repulsion respectively, and the last term corresponds to 
the bare ferromagnetic {Jh > 0) coupling according to 
the first Hund's rule. For simplicity we have omitted spin 
anisotropics arising from spin orbit coupling [I9l |43l [45] . 
The second term of Eq. ([2| corresponds to two non- 
interacting leads [u = L (left) or R (right)]: 

-^leads ~ ^ ^ukc]y]^(jCukai (4) 
vka 

and the last term is the mixing (tunneling) between dot 
and leads 

^ '^tuiclj.^di^^'R.c., (5) 

iuka 



-fornix — 



where we assume for simplicity constant tunneling 
amplitudes t^i near the Fermi energy, and the same density 
of states in each lead. 

3.2. Non- equilibrium diagrammatic calculations 

A disadvantage of the generic model ([2| is that it has many 
parameters, and shows a great deal of physical regimes, as 
discussed in the introduction (see for instance Ref. 08] 
for a detailed analysis in some other interesting range of 
parameters). The singlet-triplet quantum phase transition 
(QPT) takes place inside the N = 2 diamond, but we have 



seen in Section 2.2 that charge fluctuations towards the 
closest Coulomb diamond {N = 1 or A/" = 3) are important 
for the dispersion of the magnetic excitations. Focusing 
on the triplet side of the transition from now on, we will 
thus retain the lowest states of the N = 2 and A^ = 3 
configurations only. In fact we shall show that the resulting 
singlet-triplet Anderson model (STAM) [26 is the minimal 
model that describes the QPT as the gate voltage increases 
from the transition point towards the triplet domain. The 
states and energies retained in the STAM are the triplet 
state for A^ = 2 composed by 

111) = ^1^4^10), = ei + 62 + Ui2 - Jh/^. (6) 
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Figure 8. Upper panel: experimental temperature- 
dependent linear conductance on the triplet side from far 
to close to the transition (bottom to top). A dispersive 
shoulder associated to the underlying singlet excitation is 
observed. Lower panel: theoretical calculation of G{T) (full 
line) as performed in Section |3.2| showing also the relative 
contributions from the singlet (dash-dot-dotted line) and 
triplet (dashed line). Parameters are D = 20, Et = 0, 
Es = 0.1, Ti = 1,T2 = 1/2. 



and the states |10) and |1 — 1) obtained by applying 
succesively the lowering operator 5^ + to it, the lowest 
singlet state for A' = 2 (assuming ei < 62) 



|00) =d{^d{JO), = 2ei + 
and the lowest doublet for A" = 3 



(7) 



-e2 + t/i+2[/i2, (8) 



and|;) = (si +52)|t). 

The problem of projecting the Hamiltonian onto 
this reduced Hilbert space takes the same form as that 
explained in detail in Section II B of Ref. [69 , making an 
electron-hole transformation hlj.^ 



a\ = -d2i, a| 



t 



i2t7 



6! = 



Cukl^ ^uki 
di^. 
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Assuming tLitR2 = ^L2^i?i, which corresponds to the 
situation of only one screening channel discussed in Section 

1/2 



2.3 



and defining h, 

with/it^ = E „ 

of Eq. ([2| with now 



k ^ika/ > 



^Nk, the Hamiltonian takes the form 



Fdot = Es\00){00\+Et J2 \^M){lM\+EdY,\^)H 



M 



Hi 



leads 



vk(7 




K(4U)-/i||t)) (oo| 

Vt[h\\i) + hl\t)) {io| 
\/2T4(/i||t){ll| + 4U)(l-l|) 



(10) 



(11) b 10 



where Vs = [ti^ + ^L]^^^ and Vt = [{i 



Ll 



)/2] 



H.c, 



1/2 



We will show that our main experimental findings can 
be naturally explained by out-of-equlibrium diagrammatic 
calculations using the non-crossing approximation (NCA) 
on the STAM [Eqs. ([2|, (|9|, ([^ and ([iT])], following 
our recent theoretical work [68[ |69] (some aspects of the 
equilibrium situation are also reported in NRG simulations 
in Refs. [28, 74, 31, 32 and in Sec. 3.3 below, as well as in 
the related study by Logan et al. [38 ). Let us now review 
our theoretical results, skipping the mathematical details, 
and try to connect them to the experimental discussion of 
section [21 

Recent theoretical work by two of us [68l |69] 
has confirmed one striking experimental observation 
pertaining to the finite bias conductance, namely the 
emergence, as a function of temperature, of a three peaks 
structure on the triplet side of the QPT, see Fig. [9] 
comparing the data and the numerics (note that the data 
correspond here to a gate voltage taken on the triplet 
side of the transition in the lower panel of Fig. [T]). In 
this theoretical calculations, as in those discussed below, 
a symmetric coupling of the dot with both left and 
right leads has been assumed for simplicity. Allowing an 
asymmetric coupling manages to improve the comparision 
with experiment (see Fig. 11 of Ref. [69]). Studying 
now the temperature-dependent linear conductance, a 
shoulder-like feature is obtained from the calculation 



(lower panel of Fig. 
observation of Fig. 



), quite similar to the experimental 
(upper panel). This is explained 
by the admixture of the triplet and singlet excitations 
at high temperatures. For lower temperature than the 
triplet-singlet splitting, the conductance continues to 
increase, due to the triplet Kondo effect. We note that 
underscreening is not quantitatively captured by the NCA, 
which misses the logarithmic saturation towards the zero 
temperature conductance (not shown). 

These previous calculations have been performed for 
equal hybridization strength Vg = Vt, and in this case 
the distance to the QPT is simply controled by the 
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Figure 9. Differential conductance dl/dV as a function of 
voltage on the triplet side of the singlet-triplet transition for 
different temperatures. The experimental data [34| is shown 
on panel a, and the calculations 68 on panel b are made for 
the same parameter as Fig. [s] 



parameter 5 = Et — Eg. Furthermore, at the quantum 
critical point {5 = 0), the model is exactly solvable 
[26l |69] allowing to test the diagrammatic calculations 
[W. However, experimentally the bare splitting 6 is 
approximately constant as the gate voltage Vg is varied, 
while the renormalized singlet-triplet splitting is gate- 
dependent due to the orbital asymmetry, so that the triplet 
is lowered in energy with increasing Vg. In our calculations 
for positive 6 the QPT transition can take place as a 
function of increasing Vg if Vt > Vs (co rresponding to 



2.2 



Actually, the 



tLi > V^tL2) as explained in Section 
STAM for (5 > and Vt > Vs has also been proposed for 
Tm impurities and studied with NRG fifteen years ago [26] . 
As an example, defining 

F, = 2^^|t,,p(5(a; + e,fc), (12) 

ku 

A =Ed-{Es^3Et)/4, (13) 

with the Fermi level as the origin of energies, for T2 = 
ri/4, S = 0.345ri and a band width D = lOOFi, the 
quantum critical point (QCP) determined by NRG is 
at A = Fi. For smaller (larger) values of A, the ground 
state of the system is a doublet (singlet) and the system 
is in what was denoted as the triplet (singlet) side of the 
transition. 

To extend these calculation to the non-equilibrium 
situation and to compare with the experimental non-linear 
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conductance (upper panel of Fig. [6|, we have used the 
NCA to calculate the differential conductance G{V) for a 
new set of parameters with ^ > and Vt > Vs^ in which 
the transition is driven by gate voltage, thus shifting the 
singlet-triplet splitting as in the experiment. The results 
are shown in the lower panel of Fig. [6j The quantum 
critical point is at Ed = 4.3. For Ed < E% GiV) 

shows a Kondo peak at F = due to the partial screening 
of the spin = 1 localized states. As E^ — Ed increases 
into the triplet region a mild peak appears at finite bias 
and displaces to larger values of which corresponds to 
finite energy singlet excitations [68l [69]. For Ed > E^, 
the ground state is a singlet and the conductivity drops 
at y = to a value which depends on the occupation 
and is determined at T = by the Friedel-Luttinger 
sum rule [38l |69l. As Ed — E^ increases further in the 
singlet region, a sharp Kondo-enhanced [39 finite-bias 
conductance peak develops due to excitations of triplet 
character [69] . This feature becomes sharper and displaces 
to larger V far from the critical point, in agreement with 
the experimental data of Fig. |6] (upper panel). Overall, 
the qualitative agreement between theory and experiment 
is good. There is a discrepancy in the evolution of the value 
of G{V) at its maximum, which increases sharply as the 
experimental system is driven deep inside the singlet side 
of the transition, while it decreases slowly within our NCA 
calculation. Because our truncated STAM Hamiltonian 
considers only the N = 2 and = 3 states, and neglects 
the charge N = 1 configuration, valence fiuctuations are 
artificially suppressed in the singlet phase. Hence, the 
Kondo scale associated to the enhancement of the out- 
of-equilibrium cotunneling line is underestimated in the 
calculation, so that the finite-bias peak cannot sharpen 
by approaching the diamond edge. Inclusion of the charge 
= 1 states is infact be required for a better modelization 
of the data, as we will now see. 

3.3. Numerical Renormalization Group analysis 

We are coming back now to the complete model ^ 
which contains all possible charge states, from zero to four 
electrons, in the two-orbital Anderson model. The goal 
here is to present a more global view of the transport 
properties of such quantum dot, extending thus our 
previous theoretical analysis to the neighboring diamonds, 
as well as giving a confirmation of the results from the 
diagrammatic method within the singlet-triplet sector of 
the two-electron configuration. The technique which we 
use from now on is the Numerical Renormalization Group 
(NRG), see Refs. [151 [16] foi" a review. Interesting studies 
of the properties of model ([2| within the NRG, as well as 
from perturbative arguments, appeared in Refs. [28|l7H[38l 
[3T| [^ but within a parameter range that does not fully 
correspond to the present experimental conditions. We 
thus wish here to present an investigation of transport in 
the situation of asymmetrically coupled orbitals, leading to 



the gate-controlled singlet-triplet transition [3TJ [38] , which 
was not completed in full detail in these previous works. 

It is important to make a few words on methodology 
here. At present, a fully non-equilibrium generalization 
of the NRG is still lacking, despite recent advances |75], 
so that we will consider a near-equilibrium situation. 
This approximation is relatively well-obeyed, owing to 
the non-unitary conductance of the spin S = 1/2 Kondo 
ridge (see Fig. [T]), that result from a generic asymmetric 
coupling of the molecule to the source and drain. We 
also assume that a single channel couples to the impurity 
states, which seems also the case in the experiment, 
otherwise the underscreened Kondo effect would be spoiled 
at low temperatures. This situation can be achieved for 
instance by considering the following constraint among 
the tunneling matrix element: tLi = cos{(j))ti and tRi = 
sin(^)ti, with (j) an arbitrary angle that encodes the 
maximum value Go of the conductance in the device. 
Indeed, from the hybridization to each lead = 7rpo(^ii + 
^1^2) = ^Poit'ni + t'R2) with po the density of 

states in the leads, one gets Go = (2e^//i)4ri,ri^/(ri, + 
Tr)'^ = (2e^//i) sin^(20). One can then perform a 
linear combination of the left and right lead electrons 
to end up with an effective single-channel two-impurity 
Anderson model with tunneling elements ti and ^2 onto 
each orbital [28l [32l [38] . One can compute then the total 
spectral density (or T-matrix) within the standard NRG 
at equilibrium [28l [38] 

pH = -- Y1 r ^IV ^^G.j(^) (14) 

where the hybridization matrix reads Tij = TTpotitj^ and 
the retarded Green's function Gij = ((d- ; dj)) is computed 
on the real frequency axis. This spectral density p{uj) 
requires very high resolution data in order to describe 
reliably both the diamond edges (linked to the Hubbard 
satellites) and the fixed-charge cotunneling features within 
the Coulomb diamonds. In order to obtain the latter, 
the extensive use of an optimized broadening method, the 
so-called "b-trick" [76], turns out to be crucial (technical 
detail will appear elsewhere [70]). One finally uses the 
Landauer-Meir-Wingreen formula [TT] [38] in order to 
compute the finite-bias conductance: 



+00 



G{T, Krf) = Go / duj 7r(rii + r22)p(w) 



(15) 



^(u; + Vsd/2) + ^{oj-Vsd/2) 



where n{uj) = l/(exp(a;/T) + 1) is the Fermi function. 

Because of the near-equilibrium approximation for the 
T-matrix, we cannot pretend to a full quantitative agree- 
ment with the experimental data, as was however achieved 
for the Kondo related features at zero bias [18]. However, 
the precise determination of the many parameters in the 
two-orbital Anderson model ([2| remains a very challenging 
task, so that we can only hope to capture the qualitative 
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Figure 10. Conductance obtained by the NRG for the 
two-orbital Anderson model ([2j at zero temperature, with 
parameters Un = U22 = 0.5f/i2 = 1, Jh = 3, ti = 0.25, 
t2 = 0.125, 62-61 = 0.742 (all in units of the half-bandwidth 
of the leads), giving rise to a complex Coulomb diamonds 
diagram. One readily identifies the usual charging pattern 
as a function of gate voltage Vg and source-drain bias Vsd, 
showing that the dot can accomodate from zero electron 
(for the most negative Vg) up to four electrons. However, 
more complex structures are also noted, such as sequencial 
tunneling lines that run parallel to the N = and N = 4 
diamond edges as well as Kondo features within the charge 
A/" = 1, 2, 3 diamonds. A close-up within the middle diamond 
containing two electrons is given in Fig. |11| 



details of the experiment. An important point that we wish 
also to stress is that the observed experimental features 
possess remarkable universal features, so that the exact 
determination of the microscopic couplings in the model is 
not too important. We have already emphasized here and 
in Sec. 2.3 that spin S = 1/2 and S = 1 Kondo anomalies 



provide universal scaling function for the temperature de- 
pendent conductance, that are faithfully captured by the 
NRG calculations [18 . We will further demonstrate here 
that the gate-dependence of the singlet to triplet cotun- 
neling excitations is also a very robust property in the 
model as soon at the tunnel couplings tj^i and tjy2 
between the leads {u = L^R) and each of the orbit als are 
made asymmetric, in agreement with general experimental 
observations of this phenomenon [65 , 66", IS^jlGT. 

Let us first discuss the zero-temperature global 
transport "phase diagram" of the two impurity Anderson 



model, see Fig. [TO) Several Coulomb diamonds, showing a 
charging pattern from zero to four electrons as a function 
of gate voltage Vg^ are obtained as expected. Sequential 
tunneling lines [3 , running parallel to the edges of the 
diamonds with = 1 and A' = 4 electrons are also 
seen, and can be attributed to the level splitting |e2 — ei| 
between the two orbital states. Although these latter 
features could be anticipated from general grounds, their 
observation within an NRG calculation was not reported 
to our knowledge in previous studies [28l [32l |74] . One can 
also barely guess that this line moves within the A' = 1 
diamond as a smoother cotunneling line, an effect that 
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Figure 11. Singlet-triplet unscreening quantum phase 
transition obtained from the NRG calculations at zero 
temperature (with the same parameters as in Fig. |10| >. 
Finite bias cotunneling features, with sharp Kondo-like 
enhancement, are obtained thanks to optimized broadening 
methods [76]. The zero-bias Kondo anomaly for gate voltages 
Vg > 0.25 is associated to the underscreening of the triplet 
state, while a sharp gap occurs on the singlet side of the 
transition for Vg < 0.25. 

is also seen in the experimental data in Fig. [l] Several 
Kondo anomalies are also observed at low bias: i) in the 
charge A" = 1 diamond, the well-studied S = 1/2 Kondo 
effect occurs, with full screening of the impurity spin; ii) 
in the charge N = 2 diamond, a singlet-triplet transition 
can be already guessed, see also the close-up in Fig. [iT 
and further discussion below; iii) in the charge A" = 3 
diamond, a complex Kondo state occurs, as previously 
examined in Refs. [ST, 'SS^. In this last case, a spin 
S = 1/2 state is realized by totally filling the lowest 
orbital with two electrons and adding an unpaired electron 
on the upper orbital. Nevertheless the resulting Kondo 
screening process turns out to be subtle at the boundary 
between the N = 2 and A" = 3 diamonds. This is 
because full screening of the S = 1/2 moment provides 
a zero entropy ground state within the A^ = 3 diamond, 
while the S = 1 underscreened Kondo effect within the 
N = 2 diamond has a remanent log (2) entropy, leading 
to a second quantum critical point precisely at the edge 
between both diamonds, besides the critical point located 
within the N = 2 diamond, that we discuss now. 

We finally focus the discussion on the charge A' = 
2 configuration, for which the fine structures are given 
The singlet-triplet unscreening transition is 



in Fig. 11 



first evidenced by the obvious dichotomy in the low bias 
transport properties: a low conductance state in the singlet 
sector (gate voltage Vg < 0.25) is replaced by a highly 
conducting spin S = 1 Kondo anomaly in the triplet 
sector. Concomitantly, a clear crossing of excited states 
occurs precisely at the quantum critical point. These 
cotunneling lines show sharp Kondo-like enhancement, as 
discussed previously with respect to both the experimental 
data (Sec.|2| and the diagrammatic calculations (Sec. 3.2). 
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Clearly, the inclusion of the charge N = 1 state in the NRG 
simulations allows us to recover the correct enhancement 
of the satellites on the singlet side as the gate voltage 
is decreased towards the diamond edge. This property 
is also crucial to obtain theoretically the bending of the 
cotunneling line near the edge, that is also clearly seen 
in the experiment (see Fig. [4|. On a technical level, we 
emphasize that such numerical plot requires to perform 
the broadening of the NRG raw data with the so-called 
b-trick [73, otherwise the cotunneling lines are barely 
resolved. On a more physical ground, we have checked 
by varrying the microscopic parameters of the model that 
the gate-induced singlet-triplet transition remains robust, 
as long as the bare singlet-triplet gap does not overcome 
the energy gain associated to the asymmetric tunneling 
probabilities (in agreement with the general discussion in 
Sec. 2.2). This observation highlights the singlet-triplet 



transition as an almost generic feature of two-electron 
quantum dots. 

4. Review of various scenarios and possible 
instabilities in a two-electron quantum dot 

We would like to finish by a general discussion of 
the possible mechanisms that can drive instabilities in 
quantum dots with two electrons, in cases where only 
unique singlet and triplet states are involved. Several 
scenarios can be distinguished, whether a single or two 
screening channels are active, and depending if a Zeeman 
effect splits the triplet levels or not. The additional role of 
magnetic anisotropics due to spin-orbit coupling will also 
be considered. 

4-1- Singlet-triplet ''transition^': single channel case 

We will start by discussing the situation of the so-called 
singlet /triplet "transition", which is rather a magnetic- 
field induced Kondo effect due to the extra degeneracy 
achieved at the crossing of singlet and triplet states. 
Several proposals for this effect have been made, and we 
present first the simplest one, which involves a Zeeman- 
split triplet, in the single channel case [20, 21 . The role 
of the magnetic field here is to create an effective doublet 
by bringing the lowest triplet state in degeneracy with the 



singlet state (see panel a in Fig. 12). This effective doublet 
is then screened by the conduction electrons in a standard 
Kondo process, giving rise to a single-channel Kondo scale 
Tx, while deviations from the crossing point (a role played 
by the singlet-triplet splitting l^^^ — Et,-i\) act as an 
effective magnetic field (see panel b in Fig. [l2|. As this 
energy scale becomes larger than T^, the Kondo resonance 
is killed, and the linear conductance decreases (see panel f 



in Fig. 12). Clearly, there is no quantum phase transition 



in this problem, but rather smooth crossovers when the 
two relevant energy scales match together (see panel c in 



singlet (zero entropy state) throughout the whole phase 
diagram. This of course is true including the degeneracy 
point, because Kondo screening quenches the entropy of 



the doublet (see panel d in Fig. 12). 



We finally note that this effect has been studied 
experimentally in carbon nanotube quantum dots [20 , and 
should be observable typically in molecular devices as well, 
due to the higher g-factor compared to semiconducting 
systems. Clearly our device [34 , taken on the singlet 
side and with a strong magnetic field, satisfies most of 
the conditions for this scenario (see the similarity of 
panel b in Fig. [2] and panel e in Fig. 12). However, 
the Kondo enhanced crossing of the singlet and lowest 
triplet predicted above [21] is clearly missing in our 
experiment, because the two orbit als seem to be quite 
asymmetrically coupled to the leads, so that the resulting 
Kondo temperature is extremely low. 

4-2. Singlet-triplet ''transition": two-channel case 

An alternative scenario for the singlet-triplet transition, 
which involves different degrees of freedom (namely a 
singlet and a three- fold degenerate triplet) has also been 
considered [23l |24l [25], following experiments in vertical 
quantum dots [22]. In that case, orbital quantum numbers 
are conserved during the tunneling processes, so that two 
screening channel are active. The control parameter here 
is a perpendicular magnetic field, which by orbital effects 



tunes the singlet-triplet gap (see panel a in Fig. 13). Due 
to the presence of two screening channels, the triplet phase 
is fully screened (this remains true at the degeneracy point, 
where a SU(2) xSU(2) state is generated), so that the 
ground state entropy remains zero at all values of the 
magnetic field, and again crossovers rather than quantum 
phase transition are expected (see panels c and d in 
Fig. 13). Because the triplet phase involves a Kondo 
process, the Kondo temperature is enhanced compared 
to the singlet phase (see panel b in Fig. ulsj). The 



low-temperature conductance thus shows an asymmetric 
maximum as a function of magnetic field (see panel 
f in Fig. 13), which distinguishes this case from the 
single channel situation discussed previously. We note 
that magnetic transitions appear also differently (compare 
panel e in Fig. 12 and Fig. 13). 



4.3. Singlet-triplet unscreening quantum phase transition 

The situation that we have examined in our experi- 
ment [34 is markedly different from the two scenarios dis- 
cussed before: the gate-tuned singlet-triplet transition in- 
volves the crossing of a singlet with a three-fold degenerate 



Fig. 12). This is because the ground state remains a spin 



triplet (see panel a in Fig. 14) in the presence of a single 
screening channel. The physics is then more exotic on sev- 
eral accounts. First, the triplet phase displays an under- 
screened Kondo effect [17 , so that a zero-bias peak coexists 
with the finite bias singlet excitations (see lower panel in 
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Figure 12. Singlet-triplet "transition" in the single-channel case: Zeeman induced Kondo effect by crossing of singlet and lowest 
triplet state. 
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Figure 13. Singlet-triplet "transition" in the two-channel case: magnetic-field induced Kondo effect by crossing of singlet and 
threefold degenerate triplet states. 
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Figure 14. Singlet-triplet quantum phase transition: gate-induced crossing of singlet and three-fold degenerate triplet with a single 
screening channel 



Fig. [T] and panel e in Fig. 14). We have later on confirmed manifested by a step like behavior of the linear conduc- 



in detail this underscreened state, see discusion in Sec. |2.3 
and Ref. [18] . Because of the partial screening of the 5* 



spin, the entropy of the triplet phase saturates to log (2) at 
low temperature, so that there is an entropy change to zero 



at the singlet-triplet crossing (see panel d in Fig. 14). This 
corresponds to a zero temperature singlet-triplet quantum 
phase transition [|Mll23llSl2S[30l|3ll[3a|33], which is 



tance as a function of gate voltage (see panel f in Fig. 14). 
The situation on the singlet side close to the quantum crit- 
ical point is also intriguing, because singlet and triplet 
states are released into two independent spin S = 1/2, 
that are screened in a two-stage Kondo process. The re- 
sulting energy scales and the associated scenario are given 



on panel b and c in Fig. 14 respectively. One can also un- 
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derstand the singlet-triplet QPT by an unscreening tran- 
sition, where the underscreened Kondo phase corresponds 
to the ferromagnetic side of an effective S = 1/2 Kondo 
model, while the singlet binding is associated to a Kondo 
screening process (antiferromagnetic side in the effective 
model) [28l[3l], see the schematic flow diagram in Fig. 17 



Roch et al. 
Parks et al. 



4-4- Zeeman splitting of the underscreened Kondo 
anomaly 

We comment here on the effect of a magnetic field onto 
the underscreened S = 1 Kondo state. The standard 
(fully screened) S = 1/2 Kondo state, governed by Fermi 
liquid theory, is known to present a threshold in magnetic 
field (of the order of the Kondo temperature) at which the 
finite bias conductance exhibits a splitting [6 . In contrast, 
we have found both experimentally and theoretically [18 
that the underscreened state has a dramatic sensitivity 
to magnetic field, due to the finite log(2) entropy in the 
ground state, see Fig. [15] In that case, the non-linear 
conductance splits (at zero temperature) for arbitrarily 
small values of the magnetic field. 

4-5. Magnetic anisotropics driven quantum phase 
transition 

We end up by a discussion of possible relevant perturba- 
tions to the physics discussed above, that are related to 
spin orbit interactions. Recent experiments [41][42[|19 and 
theoretical works 031 IH |45l [HI |M1 Hi] have demonstrated 
that magnetic anisotropics can strongly affect the nature 
of the magnetic states in a two-electron quantum dot. Let 
us consider the addition of a magnetic anisotropy term 
D^{Sz)'^j which splits the triplet into unpolarized |T, O) 



and polarized |T, ±l) states (see panel a in Fig. 16). The 
experiment by Parks et al. [19 corresponds precisely to the 
situation of > 0, where a spin = 1 Kondo anomaly 
is progressively destroyed by lowering in energy the unpo- 
larized state, so that the conductance is suppressed (see 
panel f in Fig. 16). This is experimentally realized by 



tuning the magnetic anisotropy with a stretching of their 
molecule, and this realizes the schematic flow diagram in 
Fig. [T7| In contrast to our experiment where the under- 



screened Kondo was associated to a whole phase of the 
flow diagram, the Parks experiment realizes the under- 
screened Kondo as a quantum critical point associated to 
a Kosterlitz-Thouless QPT. However, the doublet phase 
corresponding to < was not accessible in this exper- 
iment, so that the QPT was not fully examined. As was 
noted previously [43l [45] , the magnetic anisotropy term is 
less relevant than the Zeeman splitting discussed in the 
previous paragraph and in our work [18], leading to a 
threshold behavior at "^k^^ splitting of the 

non-linear conductance (very similar to the Zeeman effect 
in S = 1/2 quantum dots). This effect is quite clearly 
observed in the data of Parks et al. ^19j. 



Ferromagnetic 



Antiferromagnetic 




Figure 17. Flow diagram associated to our experiment 
Refs. [181 134] and the experiment by Parks et al. 1191. 



We finally note that some of the data with 
semiconducting quantum dots with two electrons [36| |73] 
are likely to display some degree of spin-orbit effects, 
although this is not yet fully clarified (see Ref. [45 for 
a recent discussion). 

5. Conclusion 

We have shown in this review that several striking phenom- 
ena in two-electron quantum dots (gate-dependent Hund's 
rule, singlet-triplet unscreening quantum phase transition 
and underscreened Kondo effect) can be explained in a co- 
herent picture at the light of the experimental data ob- 
tained from a molecular transistor. Extensive Numeri- 
cal Renormalization Group simulations as well as out-of- 
equilibrium diagrammatic calculations were used to con- 
firm quasi-quantitatively these observations. We argued 
that many of the observed features can be considered as 
generic for molecular devices, in the absence of strong spin- 
orbit effects. One crucial and common ingredient is the 
presence of a single active screening channel in the acces- 
sible temperature range, which is granted by the asym- 
metric hybridization between orbital and leads. This is 
expected to be satisfied in most devices, and indeed gate- 
dependence of the magnetic excitations has been reported 
in all kinds of quantum dot systems: semi-conducting de- 
vices , carbon nanotubes [63j [65l |64] , several kinds of 
molecules [34l [67] . The singlet-triplet quantum phase tran- 
sition [3T imposes however a further experimental require- 
ment, namely a small bare singlet-triplet splitting, in order 
to allow the gating effect to overcome it. Recent theoret- 
ical work [38 has also confirmed a greater generality of 
this physics than discussed here. We thus conclude that 
the unscreening quantum phase transition and the under- 
screened Kondo effect constitute quite general paradigms 
to be taken for two-electron quantum dots. 
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Figure 15. Magnetic field splitting of the underscreened S = 1 Kondo state. 
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Figure 16. Spin-orbit quantum phase transition: underscreened Kondo anomaly in the presence of a tunable magnetic anisotropy 
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